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I. INTRODUCTION 

Sanece werd War ii, combat modeling, Simulation and 
analysis have been the subjects of considerable research. 
The okjectives of this research are to support defense deci- 
sion making and doctrinal developments during peace and war 
time. During peacetime defense-planners are primarily 
concerned with weapon procurement, development, acquisition, 
organisation and structuring. During war time it is 
kelieved that a better understanding of the quantitative 
aspects cf attrition can help commanders make better command 
and control decisions. 

Cembat processes involve complicated interactions 
ketween oppcsing forces. These interactions are often influ- 
enced by many external factors such aS environment, troop 
Gudlity and™ tactics. There are different types of comkat 
models such aS war games, simulations and analytical models. 
A fundamental recuirement for a good model is that it must 
ke of a fairly high degree of operational realism, since 
otherwise they would not be credible to military planners. 
Cn the other hand, excessively complicated models can make 
the mathematics too difficult tc handle. 

In this thesis, a generalised Lanchester [Ref. 1] model 
which contains area-fire, aimed-fire, self-attrition and 
replenishment coefficients is used. It consists of a system 
of 2N bilinear equations and belongs to the general category 
of analytical models. The model is rich enough to treat 
modern combkined-arms operations involving heterogeneous 
forces. It is also fossible to extend the model to analyse 
operations cn two or more fronts. 

Ameng the many important issues that could be analysed 


using this model, the problem of optimum force distribution 


had Ereen studied by tkozencraft and Moose (1983). In their 
Faper {Ref. 2], an objective function was chosen as’ the 
difference of the aggregate attrition rates. It was shown 
that the optimizaticn problem is mathematically equivalent 
to a matrix game. Hence, the model has a saddle-point solu- 
tion with corresponding optimum force distribution vectors 
x and y for Blue and Crange forces respectively. 

In addition, the neighborhood stability of the model at 
the operating point (x” and y*) was also investigated. By 
defining two parameters, Kil and K2 which are obtained by 
considering small perturbations around the operating points, 
a great deal could tre learned about stability. 

Motivated by these results, much of the work done during 
the initial part of this thesis was directed at studying the 
effect of stability cn battle outcome. The ultimate guestion 
is, how do we exploit the knowledge of stability of an oper- 
ating point to influence battle outcome? Before this ques- 
tion can te answered, it appears that there is 4a need for a 
ketter understanding of the equilibrium points. Chapfter III 


is devoted to finding and understanding the equilikrium 


solutions and their stability behavior. Like many other 
nonlinear system of equations, the Lanchester's model 
adopted here has multiple equilibria. Stability analysis 


[Ref. 3] of a non-linear system is usually done Ly methods 
which do nct reguire prior knowledge of the equilibrium 
solutions. One example of such a method is the Liapfpunov 
method [Ref. 4]. If, by some realizable means, the equilib- 
Yium soluticns can be found explicitly then there is no need 
to rely on these indirect methods which are often difficult 
to implerent. 

One of the reasons for resorting to the Liapunov metkod 
is the difficulty in obtaining equilibrium solutions of a 
non-linear system. Many numerical methods are unsuitable for 


reasons such as difficulty in ,obtaining,..good igi 
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guesses, non-convergence, Po conditlLoningsand So forth. 
Fomntumately, a powerfui numerical technigue called the 
Continuation Method can be applied for our purpose. This 
method not cnly finds all the solutions (i.e. it is exhaus- 
tive), it dces not even require initial guesses. 

In order to gain a firm grasp on the dynamics of the 
system surrounding the equilibria, it is helpful te temfo- 
Larily fecus attenticn on the homogeneous (1*1) systen. In 
Spite or its simplicity, the 1*1 system is not devoid of the 
essential characteristics of the N*N system. In fact, the 
1*1 model is sufficiently sophisticated for certain analyses 
in which the opposing forces can be assumed to’ te homoge- 
neous. AS we proceed through Chapter IV, it will beccne 
clear that much insight into the stability and systen 
dynamics could be gained by merely considering the 1*1 
system. Part of the chapter is devoted to the derivations 
and interpretations cf the relations between system asymp- 
totes, locations of equilibrium points and stability. The 
dynamics of the system are studied using the idea of fhase 
trajectcries. These trajectories represent changes cf force 
levels with time and they will be shown to depend not only 
on the stabilities of equilibrium points but also on the 
domains of attraction. 

Chapter V concentrates on Eattle outcome which is cne of 
the main issues facing a commander. It encompasses many 
issues such as, (1) Who will win and by what margin? (2) 
What is the length of battle? (3) How do initial deployments 
affect battle outcome? (4) Which parameters affect battle 
cutcome most? But we will cnly address the two following 
subjects ; 


(a) The effect cf stability on battle outcome; 
(b) The effect cf varying X and Y, the initial force 


levels. 
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The ~tasic approach is to define a nultistage battle with 
a predetermined condition for termination. The resultant 
payoff matrix can then be used to obtain the optimum set of 
mixed strategies. An example, which employs KOREAN WAR data, 
1s presented for the purpose of illustrations and 
discussions. 

The essence of tle findings are: 

1. Unstable operating conditions can be exploited to 
influence battle outcome, especially when total war 
rescurces are large. The effect on battle outcome is 
mcre pronounced for highly unstable warfare; 

2. Initial force deployment can be optimizei in accor- 
dance with a set of mixed strategies. 

We conclude this introduction by stating two of the 
outstanding issues. The first question is the extent to 
which cne can replace the N*N problem by the 1%*1 problen. 
The motivation to find an equivalent 1*1 systen stems fron 
(1) our better understanding of the 1*1 system, (2) ease of 
presentirg and visualizing two-dimensional pictures, and (3) 
Savings in computaticnal effort. 

The second question concerns replenishment rates. Dn 
this thesis, the replenishment terms used in the model have 
keen constant. It is therefore reasonable to ask, how to 
modify replenishment terms to reflect a higher degree of 
operational realism? In other words, are there more suit- 


able time-dependent replenishment rates r(t)? 
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II. IANCHESTER'S EQUATION 


mE ee eee a SE 


A. BACKGROUND 


Ccembat models have been Studied as a form of decision 
aid for defense flanring. A wide variety of defense plan- 
hing problems, ranging from force structuring and wearfon 
selection to rates of deployment in battles have been anal- 
ysed wsing combat models. There are many different types of 
models. They can be loosely categorized as either war games 
Simulaticns or analytical models. Discussions on the 
Mature, advantages and shortcomings of each can be found in 
Ret. 5D). 

Cur attention will tke focused on a generalized 
Ilanchester's [Ref. 5] model, which 1s an analvtical model. 
It consists basically of a system of ordinary differential 
equations describing the Mutual interactions ‘between 
opposing combat forces. Although earlier works Li 
Lanchester's model [Ref. 6] employed only a few terms in the 
€guaticns, modern high speed computers enable more general- 


ised, realistic and responsive versions to be used. 


Consider a battlefield with opposing forces, Blue and 
Crangeé, denoted by i! and aa respectively. The 
SIoSchIpts i, Jj refer to the type of forces such as 


infantry, tanks, artillery, etc. A generalised version of 


Ilanchester's model given by 


= o— = . _ * aes « — * 6 a x ° 
aS x 35595 a5”; a 


J q 
an 2-1) 
= = = : -C.e.. = Fels TT : 
fo ¥52Pio3 5 D193 "a 
ae a: 
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1 eee 
J = ae ore 
where 
us , a = self-attritvonm cocfficicm:. 
ce ae = area-fire attrition coefficients 
r4% oF = almed-fire attrition coe! 2s 1ements 
ee = = replenishment eoeef:1¢1 mms 


is adcpted in this tkesis. 

Note that in general I #J, implying that the force 
compositions may be different for the two sides. It is also 
possikle to extend the atove formulation to a  scenerio 
involving more than cne battlefield. 

In the next two sections, the highlights of the work 
done by Wozencraft and Moose (1983) are given. The work done 
in this thesis 1s a continuation and extention Of ties 
Wolle ies The detailed derivations of the results obtained by 


them can be found in[{Ref. 2], and hence are not included 


Bo. OFTIMNOM FORCE DISTRIBOT IGN 


The guestion of optimum force distribution arises in 


combined-arms operations. The problem is fundamentally 
Chase Given aggregated forces X, gs how Should one 
distribute them among the different types x, and vs ee i= 
Tp le ccgle.  } = 91, 2 fee Since loss rate is one of the 


fundamental concepts in combat modeling, it is reasonable to 
choose this measure as a Starting point. The okjective 


function was chosen to be 


1? 


A ° ° 
M ay Gs, - Vy? - 4, i 22) (eqn 2.2) 
i J 


Morueetiuo meChOomGe Ot M. Lt was Showh that there exists 
optimum force distrikution (row and column) vectors x* and 


y 


“ such that for any cther vectors x and y 


nw * * sk A 
oe KAY (Ge 1G) ae roy) 
where 
* mo 6s 
Mo o= x Ay 
ani ome desormined by attrition coefficients and 


the aggregate force levels X and Y 


The resemblance of this result to the Minimax theorem 
{Ref. 7} in matrix games 1S very striking. Indeed, this 
result holds precisely because M can be written in a form 
mathematically eguivalent tc a matrix game. Conseguently, 
it iS net surprising that one can solve for the optinun 
vector x* and y* by means of a Linear Program. An interac- 


tive program to solve a 2*2 program is given in Appendix A. 


C. NEIGEBORHOOD STAEILITY 


Equilibrium conditions can be achieved if the reprlenish- 


rent rates are chosen to make 


~ 
1" 
Ke 
i 
Oo 


ao 
\a Sb 
NON 
Va 
an 
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at Xx = x* and y = (yee Following the usual approach in the 


analysis of nonlinear system stability, equation 2.1 can 


then be transformed into a system of linear equations. 


Sx 6x A B 
<a i -C s C = A ~~ 
Oy Oy be 


C is called the conflict matrix and its elements are dememe 
Mined by the attriticn coefficients and the optimum vectors 
xoeera Ciay ane For the system of eguations 2.1, A and ¢ are 
diagonal matrices. It was shown that two parameters kj, k, 
partially characterize the stability of the system. k, and k, 
turn out to be the column sums of the left and right side of 


the matrix 


_— = 2 ee ele 


aA 


“A a “A A 
Tenoting the elements of the submatrices A, B, C, D, by aise 
Bij. C5 and d,; respectively, k, and kK, can be written as 


ieee yoni 
l 

i = oe ye: | 

2 33 13 
ih 


independent of the columns i, j. Furthermore, it was shown 
that the following relation holds 


OX - kK, ox oro) 1 ies k,6Y 
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where 


A 
a ome | OY Sy; 
c ) 


* 5 ) 


It was found that the eguilibrium point (x*, y 1s 
Stable if k, and k, are negative. If k, and k, are positive, 
Paemecheysystem 2S *‘unstabdle'.i Furthermore, values cf k, 
and k, and hence the Stability of the operating point was 


meuma to be affected Ey the aggregate X and YX. 


IMore generally, it can be_shown that k,<A,<k,, where A, 
1s the maximum eigenvalue of -C. 
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III. MULTIDIMENSIONAL (N*N) SYSTEM 


rr ee ee = == a 


A. NATURE CF N*¥N PROELEM 


The interesting results highlighted in the last chapter 
provided motivation tc extend the body of knowledge. A study 
of the effect on stability of battle outcome seems to have 
important potentials for applications. Should a ccmmander 
strive to establish a stable operating point, and if so, 
under what conditions? Also, what is the optimum initial 
level of forces he should deploy and how many should he 
Maintain in reserve? To anSwer these guestions, more knowl- 
edge abkout the nature of these eguilibria and their 
Stability behavior is required. 

The next Section outlines the kind of problems we would 
expect to see and their potential complexity. Tt ee 


followed by. a section on finding the eyuilibrium soluticns. 


1. Existence of Kultiple Equilibria 





An N*¥N system is in equilibrium if the replenishment 


rates ©. , §S are such that there is no change in the force 


nee: 
levels (x, = aa QO). The system of equations becomes 


0 = es Xs Da - Baa Os 
j 


(eqn 3.]} 
Dry Dy | Deady tS 
ih 


i, j= Glee eee 


2 





where, for simplicity, 1 andj are each assumed tc have N 
types cf forces. 

A 2N-tuple vector, z 4 (Xx, y) which satisfies eéqua- 
tion 3.1341s an eguilibrium solution. Like many nonlinear 
Systems cf equations, equations 3.1 have more than one egui- 
Zeprium point. Geometrically, these equilibrium points are 
at the intersecticns of a set of hypersurfaces in the 
2N-dimensional space. To help in visualizing the geometry, 
we can look at an example using a 1*1 system as shown in 
figure 3.1. In this case the hypersurfaces simply reduce to 
hyperbolic curves. 


Figure 3.1 Equilibrium Points at Hyperbolic Intersections. 


The existence of multiple equilibria makes the anal- 
ysis of the N¥N problem very interesting but difficult. [In 
chapter IV, some illustrations cn how the locations of these 
equilitria affect phase trajectories will be presented. 

A few other interesting questions arise sSpontane- 


ously. For instance, how many of these equilibria are there 
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Figure 3.2 Infinite Number of Equilibria. 


in an N*¥N problem? ‘The answer to this question is not imme- 
diately obvious just by looking at equation 3.1; however, it 
emerges guite naturally when the Continuation Method is 
considered in Section IIIB. It will be seen then that an 


N*¥N system has, in general, N, equilibrium points where 


Two excepticns, or degenerate cases, have been 
observed, nanely: (1) when some or all of the hypersurfaces 
merge there are an infinite number of equilibrium points, 
(see Figure 3.2), (2) when some or all of the hypersurfaces 
intersect in such amanner that repeated equililfria are 
formed, the number of distinct equilibria is less than QW,. 
Figure 3.3 illustrates such a degeneracy. 
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Figure 3.3 Repeated Equilibria. 


foe stacey eamaeDOMains of Attraction 


SE —/> noe ee == ee 


Each equilibrium point in an N¥*N system may cr may 
not ke stable depending on whether or not its eguilifriunm 
point can be maintained. The property of neighborhood 
stability is important because it has a strong influence on 
the phase trajectories. Generally, if an operating foint is 
stable (the maximum eigenvalue is negative), then any 
perturbation away from that point results in the systen 
returning to the same point. Conversely, perturbations 
about an unstable foint results in divergence from that 
pomnt. 

The notion of domains of attraction is also critical 
when determining phase trajectcries. Any operating point 
within this domain cr region will be "attracted" toward a 
Stable eguilibrium point. Pieswenme,na cola of attraction 
is a volume in the 2N-dimensional space surrounding a stakle 
eguilibrium point. Figure 3.4 shows a typical domain in 
which seme of the trajectories are shown converging toa 
Stable eguilibrium point. 


24 


a 
y boundary 
4 
4 
7 
4 
/ 
4 
Ag 
4 a a 
? 4 : 
¢ W_lsé=OC«srajeectory 


\ X Unstable. vote 


\ 

\ 

: * Stable point 
t 


Figure 3.4 Domains of Attraction. 


Domains of attraction are separated by boundaries 
which are invariant curves in 1*1 problems ani invariant 
hypersurfaces in N*N problems. A boundary surface may be 
considered as an infinite number of invariant curves placed 
Side by side. A boundary curve is the locus of points that 
approach an unstable point from both sides. The boundary 
line can be obtained by backward integration (i.e. using 
negative time in equation 2.1) starting just on either side 
of an unstable point. The rationale behind this method is 
that to approach anunstable equilibriun, a point must 
remain exactly on tke boundary. If this is not’ the case, 


then the point will ke attracted into the domains and move 


2 


toward a stable point or infinity. By performing a tackward 
integration, we are actually retracing the path taken bya 
point which previously approached the unstable eguilikriun 
point. This method requires knowledge of the unstable equi- 
libria, but this 1s made feasible because the Continuation 


Methods can be used tc find all eguilibrium solutions. 


BE. FINDING THE EQUILIBRIUM SOLUTIONS 


To obtain a set cf equilibrium solutions, one _ has to 
solve eguation 3.1, which can be written using a more 


compact notation as 


Baez) = 0 (eqn 3.2) 


where 


F(.) represents the right-hand side of 
equat lon Seal 


N] 


= OS 


QO = zero vector 


It is well known that numerical technigues for solving 
nonlinear eguations are not always successful. Since equa- 
tion 3.2 describes a bilinear system, one Should expect to 
face similiar difficulties when attempting to solve it 
numerically. 

Most numerical methods for POCmetbmadrng generally 
require that a fairly good initial guess (Z,) be known so 


that some convergent iteration process 


ent+1 g(2,) 
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briugs the apprcximated root closer and closer to 7 the 
desired eguilibrium solution or root. in practice), “"Yire 
following difficulties are often encountered 

(1) The convergence condition of the algorithm 
must be ensured ; 

(2) Finding an initial guess that is sufficiently 
close to the correct solution is diffticuwiee 
especially for higher dimensions; 

(3) Even if a good initial guess has been obtained, 
the numerical process may still be plagued by 
ili-~conditioning, saddle points, etc.; 

(4) Not all the solutions are guaranteed to ke 
foun as. 


1. Continuation Method 


Fortunately, the above problems are avoided if a 
Numerical method called the Continuation Method [Ref. 8] is 
used. This technigue, which is sometimes called the 
Imbedding Method, has been successfully applied in Many 
fields. It introduces an artifical guide which will channel 
the iterates toward a specific solution. Such a  gilaig@amig 
principle is actually a knowledge of the existence of a 
Suitable curve connecting an initial point with the desired 
solution. | 

Continuation Method has significant advantages over 
other numerical technigues. Most importantly, a good initial 
guess if not necessary and all the solutions can be 


obtained. 
a. Basic Theory 


Given the problem F(z) = 0 to solve, the first 
step is to embed it into a homotopy or a parameterized set 


of problems, H(Z,t). The requirements on H(z,t) are: 
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Cie) s=8F> (z) 
(2 ) H(z,0) = Fo (Z) 


i enemomrgmnal problem 


i 
Cl ol 


has a trivial or easily 


computed solution 


For example, a homotopy could tLe ; 


epee chia (t-t)yr (z) , te[o,T) (eqn.3.3) 


Using the above parameterization, the simple 
problem of Fo (2) = 0 is deformed into the desired one, F, (Zz) 
= 0 - This is done by calculating the solution te the 
deformed problem at each stage of the deformation. The exis- 
tence of a continuous curve such that H (z(t), t) 1s a solu- 
meron to H{.,.-) = 90 for all te (0, 1] is assumed. 


bk. Implementation 


To actually carry out the above continuation 


process one usually differentiates H(.,.) to form 


H(Z(t),t) = 0 (eqn 3.4) 


Using equation 3.4, 2 can ke written as a function of z and 
t as given in eguaticn 3.5. The function, h(.,.) is prefer- 


ably a linear functicn that can be integrated numerically. 


2 = h(Z,t) (aqme 3-5) 


Together with the initial condition z(0) =z, , 
equation 3.5 1s actually an initial value problem which can 
be integrated numerically. The solution at t=1 is then the 


solution to the original problem F(z) = 0. 


2G 


2. Algorithm to Obtain 2*2 Equilibrium Problen 


SS a oS 


A 2*2 Lanchester problem is first formulated into a 
Continuation process. It is followed by a discussion on how 
the accuracy of the method can be improved. The last part of 
this subsection includes a note on the number of equilikriun 


Foints in an N*N problen. 
ae Formulation 


For the 242 problem, F(z) = 0 itgpexpl demu 


~ 2,00) 48) 123*8p cg), 9D), ee ee 
ho (Uy ii 2 5485524) tea: Oat ee 
~25(U34C1 4297059720) 7 5, “Sigal 
2g (U4*Cy 92447229) + Ty -dy 52) - dy,2, = 0 


The horotopy is formed by writing 


H, (2) = H, (2,9) 4 t(r,-b) 423-6) 524) = 0 
H,(z) = H,(z,0) + tt (3) Deeper eee 
Z Ee “ Ze te 4 (eqn 3.6) 
Hz (Z) - H,(2,0) + t(rz-d,,2,-d5,25) = 0 
H, (2) = H, (2,9) ti t(ry-d, 52,745,525) = 0 
where 
H(z, 0) - -2) (u, +4) 425481924) = 0 
H(z, 0) = "25 (U,+85,23485524) = 0 
Hz (2,9) -23(Uz+C) 2,705,275) 0 
H, (2,9) ~ -Z4(UgtC) 42, 4C5975) = 0 


ae 


Next, we differentiate equation 3.6 
puter t Inp@av matrix form 


With respect to t and 


AZ ws) (Ceicgn oer?) 
where 
Uy +4) 123444 224 0 2442, 7D jt Ay 924 +O ot 
Y Bogie 774 27)70*D2,t seme 
A = 
11" ee 171° 2172 : 
Cyotgtdyot = Ca9%4*Gz0t 0 Ba S1272°%2272 
= : i 
= = [24125523524] 
T,)-b,4235-by 974 
i eS Oy 2 
B= Vee Ne OP lcs 
T3745 471799) 
ee 5°22 





Equation 3.7 can now te integrated numerically using one of 


the readily available integration froutines. 
We have assumed that the trivial solution to 
H(z, 0) = 0 has been previously found. 


be. Improving Accuracy 


Numerical integration of eguation 3.7 inevitalkly 


produces scme errors at €ach iteration. Since the 
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Continuation method relies on following curves to arrive at 
the desired solution, it is esssential that each iterate 
remains close to the actual curve. It is necessary to 
include a way to correct the approximated position Ey means 
cf a corrector step. The combination of integration and 
correction is often called a "predictor-corrector step" 

This process of prediction-correction is Siew 
in Figure 3.5 where ¢ach integration error has keen exagger- 
ated for illustrative purposes. The algorithm tc ae 
presented later emplcys an IMSL routine called ZSCNT for the 
predictor step. Other forms of curve following routine can 
also ke found in the literature, and are briefly mentioned 
in [ here, oak 


Cc. Trivial Solution 


The trivial system H(z,0) = 0 was chosen to be 
H, (2,0) = -Z, (U, 44,4 2344; 924) = 0 
H,(z, 9) = -Z2,(u,+a,,2348,524) = 0 ‘ean Su 
H,(z,0) "22 (U34C) 424 4C 5,29) = Q 
H, (2,9) = “24 (Ug) 924405929) = 0 


In non-degenerate cases, there are six scluticns 
corresponding to equation 3.8. The result is derived in 
Appendix B which also deduces the number of trivial solu- 


tions for an N*N problem to be 


NN Z 
Ny S) (eqn 3.9) 


i 
1=Q0 


The method of oktaining the trivial solutions is given in 
Appendix 8B. Using a combinatorial identity, Ny can be 


written as 
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Figure 3.5 Integration Path using Predictor-Corrector. 


*") 
N. = 
k 
N 

Hach COmmnuadtlonm process starts from a trivial 
Soluticn Zo and follows a specific curve until it reaches 
the equilibrium point. Conseguently, the number of eguilib- 
rium points will alse be Ny» AS mentioned in section IIIA, 
the two exceptions are situations involving infinitely-many 


and repeated equilibria. Situations involving degeneracy 


are discussed in Appendix B. 
d. Algorithno 


(1) Singularity Treatment. iieeonteianitat10n 
Method algorithms [{Ref. 8], it is sometimes necessary to 


give special treatment to cases in which the curves being 
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followed by the integration reutine pass through a singu- 
Lat IMity < Experimentally, it had been observed that in our 
problem, the Singularity took on the form shown in Figure 
3.6. Corrective measures were necessary to ensure that upon 
crossing the singularity, the large magnitude was preserved 
but the sign was changed; otherwise the curve might terri- 
nate at an eguilibriug point which was not the intended one. 
In the algeria, the presence of the 
Singularity is detected by monitoring the rate of change of 
the individual component z,. Once identified, this fast- 
changing and large-magnitude component (Z,) is monitored at 
each step t where 0 = ty <.-.<.-- t, < tual Kee eS tong ‘le 
When zZ, is found not to cross the singularity and end up at 
approximately -z,, the algorithm attempts to Correct spa 
irregularity by artificially making Zn = WZ before the next 
predictor step commences. 
(2) Flowchart. The flowchart for the algo- 
rithm is given in figure 3.7. Only the major steps have been 


shown. The program listing iS given in Appendix C. 
3. Example and Fesults 


Consider as an example a 2*2 problem with the 


following attrition coefficients 


eee 022 
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Figure 3.6 Curve Passing through Singularity. 


The trivial solutions are first computed and serve 
as one of the inputs to the program. The program obtains the 
values of z(t) and plots each component ( 2; (t) ) versus t. 
In Figure 3.8, the plcts for t close to zero show one set 
mamourves LOL Z.{t) Starting from their respective trivial 
soluticns. PiemGtimicommOon Z(t) | versus t€ for all the six 
sets cf equilibrium solutions are shown in Figure 3.9. 

A few interesting features of the continuation 
process are worth noting. For example 

e Each trivial solution leads to different eyuilib- 
riun solution andthe integration path is different 


for each component. 


e All the curves are smooth ; one of the four curves 
May pass through a singularity. ( see Figure 3.9 (c) 
and (d) ). 

Table I Summarizes the computed eguilibrium  solu- 


tions. They are tabulated in the same order as the plots in 
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Pigure 3.7 Flowchart for Algorithm toe find 2°2 Eau ipa 
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Figure 3.8 z(t) Versus t for Values of t Close to Zero. 


Figure 3.9. To estimate the accuracies of the results, we 


defined error as 


2 2 
ERROR = Fi. + FL. + FL + b 
if 3 4 


: F(z) » Z 1s the computed equilibrium 


Soutien to F(z je= 0 


DecreaSing the 1notegration step size in the 
predictor and correctcr routines may reduce the errors bya 
small amount; but tke increase in computational effort ray 
not be justifiable. Conversely, it may be desirable to cut 
down computing time. Currently, the algorithm performs one 


corrector step for ¢ach predictor step. If two or nore 
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Figure 3.9 Plots of z(t) Versus t during Continuation Process 
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TABISE.. 2 
Computed Equilibria for X% = ToC Qi 0 





Triviale soluimion Computee Equilitrium selumignesz Eiri 
on (1.7755, socom 0.2° 10 
0, 0) 1.0533, 0.3512 
(GO. a S0me ee (34.8989, -Gul coum 0.610 °° 
Owes 30) =0.1350, 20sec 
(0, 20: 09mm (0.6154, 0.3946. O.1"*10°° 
=> oc 3.0769, 0.82231) 
(-0.2, 0, (0.Weeto OT See 0.22"100" 
Oxo =16-6400, 52.98z2Gm 
(-0.083, 0, (-0.1857, 0.09325, 0.56 ssnaneee 
03, oD SU4.5297, 28 -Gaee 
(-0.42, 0.37. (-45.2284, 39.1604, 0.23%107° 
0.25, -0 ae 0:07, =0. 00S) 
predictor steps are done for each corrector step, scne 


computational effort2 can be saved. 


2Saving in computational effort will be more significant 
when solving higher dimensional systems. 
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IV. PROFERTIES OF THE 1*1 SYSTEM 

The 1*1 problem is the Simplest case in our nodel. It is 
nevertheless important for us to investigate and understand 
its properties. Despite its relative simplicity, it is by no 
means uninteresting. There exists many situations which can 
ke realistically and easily modeled by the 1*1 sySten. PONS 
example, when the opposing forces can be considered as hono- 
geneous, it 1S convenient to use the 1*1 model for analysis. 
It is also useful for the analyses at the strategic level 
when the forces and parameters can be aggregated. In many 
instances, it seems to provide insight on how to approach 
the N¥N problem, which is much more difficult to visualize. 
In fact, as the understanding of the 1*1 system increases, 
there is a strong urge to try to represent the N*N prcblenm 
by an equivalent 1*1 problem. The eguivalent representation 
is not only attractive in terms of its simplicity but also 
its econcmy in computational efforts. 

The next section will focus on the relation between 
System asymptotes and stability of the equilibria. By formu- 
lating the problem quantitatively, we are able to arrive at 
some useful properties. In Section IVB, the system dynamics 
1-e. the changes in the force levels are analysed Ly consid- 


ering the phase trajectories. 


A. SYSTEM ASYMPTOTES AND EQUILIBRIUM POINTS 


For the 1*1 rroblem, the system reduces to 


x = -x(u saclay 24) eeeceanTas coe 1), 


: ecanimear 
oa eer CK ye S - dx 
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An equilibrium condition exists if r and s are chosen such 
that x and y are both zero. In general, there will be two 
equilibrium points corresponding to two locations where the 


two hyperbolas intersect. The hyperbolas are descriked ky 


y by 

u + ay 
(eqn 4.2) 

Ss - dx 


yO eae 


From equation 4.2, one can easily deduce the fouc asymptotes 
(two vertical and two horizontal) associated with the hyper- 
bolas. Figure 4.1 shows a typical set of four asymptotes. 
They always cross in the third guadrant of the x-y plane and 
do not depend on the replenishment coefficients. The rela- 
tive displacements ketween the two horizontal (and also 
vertical) asymptotes depend only on the ratios of attrition 
coefficients and not on the coefficients themselves. et 
turns out that these properties of the system asynptotes 
help to simplify the analysis considerably. 


Considering small perturbations about an eguilifriun 


(x. .Y,) and linearizing the equations, we have 


Sx (u+ ay.) (b+ ax,) Ox 


Sy Ce a 


The characteristic polynomial is simply 


u 
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Figure 4.1 System Asymptotes. 
D(s) = Det [sI - C] (eqn 4.3) 
where 
fe=eidentity matrix 
C = the 2*2 matrix in equation 4.3 
Hone, 
D(s) = E +, (utay,)|| s ee (v+ex,)| : j(b+ax,) (d+cy,)| 


sf + | (utay,) if (vtex,) | s + |(utay,) (vtex,) 


- (b+ax,) (d+cy,)] 


The conditions for (Xo ¢ 
iee. for the roots of D(s) to be in the Left Half Plane 


VeeOmbemmdescable equilibriun, 


(LHP) are given by 


4 2 


Ge ay.) + (Ve CX) > 0 


(eqn 4.4) 
(Cu ay.) (v + cx.) = ea ax.) (cia cy 4) rane 


2. Stability and Asymptote 


= = toe 


Five different ways in which the hyperbolas can 
intercept have been identified and their stabilities 
accounted for. These five cases are shown in Figure 4.2 and 
each case will be elaborated upon subseguently. 


ae Definitions and Formulation 


One of the most intriguing facets of the 1*1 
problem is the connection between the asymptotes and the 
Stability cf the resulting equilibria. We begin the quanti- 


tative treatment by first defining the following ratios: 


A ou Ad 
ee = Nn set 
1 a 2 Z S 
Ab Ay 
uo=w# yost 
i a : 2 Cc 
The four asymptotes are x = 2 XS, ey = =e and y = 


=o If we let the first equilibrium point be (X07 
and substitute the corresponding r and s into the equation 


G.2, we have 


Ul 
SS 


Ny Ox- X04) : CS Fe oi) 7 Wy YY ey) 
(eqn 4.5) 
nN, (x-x 64) a (xy -x ee “ Wo Oy-Yey) = 
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case (a) case (b) 





case (c) case (d) 


Wns Gabe 


Ss *) stable 


~ : 


case (e) 


Mewmterally stable 


Fisoune 4-2 Types of Equilmria 
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Next, the distances rFéetween the asymptotes are defined as 


Le 


€ 
x 2 1 


1 


“y pe 


It 1s not ditimcwlt tomsee tier — ead Ey W ig 
decide where the hyperbolas intersect. For instance, when €, 
> 0 and a > 0, there may be two equilibria in the first 
guadrant? (See case (b) of Figure 4.2). In general, the 
e2/ te2 
nating y or xX from equation 4.5 and comparing coefficients 
with (y- ae eur = ie) and (x - X 3) (X = Xo) - The final 
expressions are 


second equilibrium point (X can be found Ey elini- 


€ 
= x A 
eve Ey OE : Nh a 
: (eqn 436m 
5 ae : 
Ye2 Ey (Xo . ae A 


For constant %.)4 X.54 Yo ,0. and Yoo equation 
4.6 can be written to represent two straight Jines ingieaae 


Ey, plane. The equations of these two lines are 


€ c er + Ny) 


ye 
(X21) ae 
(eqn 42% 
€ € O71 - ae 
Y =O X42 eae 
(X.5 a Wy) 


3In our context, the quadrants are defined by the asynp- 
totes and not by the x, y axes. 
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Ceo Sel SEGUE bDEia and their Stability 


To derive the different types of equilibria and 
their associated stabilities, we make a transition from the 
Wa Plate into thewe,, E€. plane. Briefly, the basic 


yy, 


approach is to £1x ore equilibrium point (Xo, 7 Y on the 


) 
ea 
first guadrant hyperbolas and consider the regions in the 


xe «Ey plane when we have the other point in various places 


of the x, y plane. The other essential step is to express 


e 


the stability criteria (equation 4.4) in terms of € = 


x @ ee 
Wye Tye ¥,¢ Yo, ° A Summary of the results which are derived 


in Appendix Dis given below ;: 


(1) When both equilibrium points are on the first 
gQuadrant hyperbolas ( case (b) in Figure 4.2 ), cne 
will be stable and the other unstable; 

(2) When one equilibrium point is on the first 
quadrant and the other on the third, both can be 
unstable or one will be stable and the other 
unstable ( case (a) in Figure 4.2 ); 

(3) When both eguilibrium points areon the third 
quadrant hyperbolas, both are unstable ( case (c) in 
Faguisee 4.2 ); 

(4) When there are infinite number of equilibria as 
my Case (@)) “In Figure 4.2, €, = 2a O and the two 
sets of hyperbolas merge. Eguilibria lying on the 
first quadrant hyperbola are neutrally stable (one 
elgenvalue equals zero) and those on the cther 
hyperbola are unstabkle; 

(5) When there are repeated eguilibria as in case 
(e) in Figure 4.2, they are neutrally stable if the 
hyperbolas touch in the first quadrant ; otherwise 
they are unstable. 
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Most of the above results are embedded within 
Figure 4.3 which is reproduced from Appendix D for conven- 
lence. Evidently, kEcth the coordinates of the equila beiga 
points (X,, Y.) and the location in the Exe €y plane deter- 
mine the stabilities. The Cay om plane has been subdivided 
into a few regicns each with daystta nct stabil 


Characteristics. 
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Figure 4.3 The €,, E, Plane. 


The case of infinitely many equilibria corre- 
xe ee plane (H = use Nn, =1,)- 4tThe 
only way for two sets of hyperbolas to merge 1s for their 


respective asymptotes to merge. This case is a degenerate 


4] 


instance of repeated eguilitria ( case (e) in Figure 4.2), 
which is shown in Appendix D to correspond to operating 
points on the line Se = E (Y + n,)/(X + by) as illustrated 
ime 5 voune 4.3. 

As a corollary, we note that there cannot be two 
Stable equilibrium pcints in the 1*1 problem. This deduction 
Can be made by referring to Figure 4.3. There iS no region 
in the €,, Ee, plane which allows for this cas2. At most, 
there can be two neutrally stable equilibria which are 
repeated. Numerous attempts have been made to obtain two 
stable equilibria in the 2*2 preblem, but in vain. Whether 
it is also true for 2%*2 or higher dimensional problems that 
Only cne equilibrium may ke stable is still a matter of 
conjecture. 

In Appendix E, the relations between the regions 
on the c,, os plane and their associated stabilities are 
verified. some representative points on the «,, Ey plane 


are chosen and their stabilities checked. 


Bo SYSTEM DYNAMICS 


The dynamics of a 1*1 system are characterised by its 
phase trajectories, which are curves on the x-y_ flane 
descriking the history of the system as the time, tee 
changes. These trajectories can be conveniently obtained by 
integrating equation 4.1 numerically. 

Needless to Say, being able to predict the trajectcries 
is important, for it means that we know how our model of a 
Fattle progresses. Cnce the factors influencing the course 
of a battle are known, appropriate command decisions can be 
introduced to ensure favorable battle outcome. In Chapter V, 
we will see how many cf£ the results obtained in this section 


can ke used to rationalize and predict battle outcome. 
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some typical trajectories corresponding to the different 
types of equilibria are described in the next subsecticn. 
Besides the stability which influences trajectories, it was 
briefly mentioned in Chapter III that domains of attraction 
also affect the trajectories. In the subsection that 
follows, we will show specific examples of the way to deter- 


Mine the domains by finding their exact boundaries. 


1. Trajectories 





Two methods cf establishing the trajectories froma 
given initial condition will be described. The brute-force 
method which has been mentioned uses numerical integration. 
The cther method which often provides better insight, is 
more graphical. The graphical method is based on a few very 
Simple rules to predict the gross behavior of a trajectory. 


Some cf these rules are listed below : 


(1) A stable point “attracts"; unstable (| pGmm: 
winewet Sa 
(2) Points an either side of a boundary move into 


their respective domains; 

(3) For large (x,y), trajectories are governed by 
the Lanchester "linear law"; 

(4) Points near the hyperbolas can be easily 


analyzed by noting the signs of x and y. 


AS an example of using the graphical method to 
determine trajectories, consider a region around an unstakle 
equilibrium point on the first quadrant hyperbola. The whole 
picture of the phase trajectories (sometimes called phase- 
plane portrait [Ref. 9] ) can be put together in a logical 
fashion by using those simple rules. Since this eguilifriun 
point is unstable, trajectories will be expectei to diverge 
ieio@) Mew amie As an unstable equilibrium point, it will havea 


boundary line passing through it. Initial conditions start 


4g 


from each side give rise to different trajectories. Next, we 
determine the signs cf x é y on both sides of each hyperbcla 
aS indicated in Figure 4.4 where only one intersecticn is 
Shown. 


X <0 
Y>0 
hyperbola 


ny 





Figure 4.4 Analytical method of predicting trajectories. 


Note how predictable these trajectories are. if, 
for some reasons, the exact trajectories are required, we 
can resort to the brute-force method. The methods are obvi- 
ously complementary in nature. The advantages of the brute- 
force method are accuracy and Simplicity. In Figure 4.5, a 
typical computer plot consisting of ten trajectories is 
Shown. The program which produces the plot is included in 
Appendix fF. 

Referring to Figure 4.5, the trajectcries cross the 


hyperbolas and move asymptotically along a common curve 
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Figure 4.5 Computer Plot of Trajectories. 


lying between the hyperbolas. This same property is exhib- 
ited by other cases. Even the special case with no hyper- 
Folic intersection has been found to behave similarly as can 
seen in Figure 4.6. 

Our ability to determine the trajectories and 
present them vividly is partly due to fact that two- 
dimensional pictures can be easily drawn and visualized. For 
dimensions higher than the third, it is impossible to visu- 
alize trajectories; however, the notion of trajectcries can 
be conceptually extended to n-dimensional space. Thus, it 
seems likely that in the higher dimensional systems, trajec- 
tories cross hypersurfaces and move along a common asymp- 
totic curve analogous to that in the 1*1 systen. Further 


Studies are required before this behavior can be confirmed. 
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hyperbola 
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trajectory 


Figure 4.6 Trajectories when Hyperbolas do not Intersect. 


2. FPoundaries of Domains of Attraction 


EE SE EE ee = det ee EE —_ sa = ee Ee ake Se ee 


In Chapter III, the idea of the domains of attrac- 
tion was briefly discussed. In an n-dimensional space, such 
a domain is a region or volume in which ali initial pcints 
come under Similiar influence. When domains exist, there 
will be boundary surfaces which can be thought cf as 
collections of invariant curves passing through unstable 
eguilibria. 

For a 1*1 problem, domains and boundaries are nct at 
all abstract. In the last subsection, they have been shown 
to affect trajectories. Recall that in Chapter III, we 
menticned asimple and yet effective way of finding the 
boundary curves and establishing the domains in the x-y 
plane. Examples on tke use of backward integration to obtain 


boundary curves are now presented. 
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a. Boundary Curve through an Unstable Point 


Starting from an unstable point, we apply small 
perturkations in both directions perpendicular to an eigen- 
vector asscciated with the most positive eigenvalue and 
integrate backward in time (in the computer program, this is 
eaSily done by employiny negative time steps for integra- 
tion). The result is a smooth, invariant curve which is 


exactly the boundary or the so-called separatrix like the 
one shown in Figure 4.7. 


Boundary 


Figure 4.7 Boundary Curve through an Unstable Point. 


To verify that the curve is indeed the boundary, 
two initial points are chosen just off the curve (e.g A, B 
in, Figure. 427) - If we forward integrate from these two 
points, they move into different domains as indicated in the 
Same diagran. Appendix G contains a Fortran program that 
does the backward integration and plots the boundary curve. 
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Boundary curves do not necessarily pass’ through 
unstable points. Backward integration methods can also be 
used if a boundary exists but there is no unstable eguilib- 
rium point to serve as the Starting point of integration. 
This is kest illustrated by considering the case of bcth 
equilibria cn the first quadrant hyperbolas. In this case, 
there is no equilibrium point in the third quadrant; never- 
theless a boundary does exist between the third-guadrant 
hypertolas. The existence of the boundary is visible by 
Simply ccnsidering tke signs of Xx and y on both sides of the 
hyperbolas. In figure 4.8, the signs of x and y and also the 


directions of some typical trajectories are depicted. 





Figure 4.8 Existence of Boundary Between Two Hyperbolas. 


To obtain the exact boundary, choose a point clcse 
to a hyperbola and cn lower part of the hyperbolas (e.g. 
point P or Q in Figure 4.8) and integrate backward. The 


result 1s a boundary curve as Shown in Figure 4.9. 
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hyperbolas a7 


Figure 4.9 Exact Boundary Curve Between Two Hyperbolas. 


4. Summary of the 1*1 Problen 


SS a ME 2. === 


We have seen the close relation between system 
asymptotes and stabilities. Through the use of newly defined 


variables c¢, ande the stability of different types of 


‘ 
equilibria has _— derived. Five cases have been identi- 
fied, and they correspond to the types of intersections on 
the x-y plane. For example, if both the equilibria are found 
ohn the third guadrant hyperbolas, then we know that tkey 
will be unstable. 

Iwo methods cf establishing the trajectories have 
been descrited in this chapter. These two methods complement 
each other and the chcice depends on our requirements. The 
dynamics of the system are characterized by the trajecto- 
ries, which as we have seen are very predictable. These 
trajectories are influenced by the stabilities of equilibria 
and domains of attraction which are separated by kEoundary 
curves. A Simple way of plotting the boundary curves has 


also been presented along with specific examples. 
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The results derived in this chapter will be applied 
in the next chapter. The knowledge of the system dynamics 
and how they are affected by stability and other parameters 
will enable us to analyze changes in force levels as the 


battle progresses. 
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In the last two chapters, emphasis has been placed on 
establishing the mathematical framework of the systen 
dynamics and stability. In this chapter, we examine some 
model operational problems that are related to stability and 
dynamic considerations. 

One of the major command decisions that has to be made 
during a build-up period of a war pertains to initial force 
commitment. A good strategy calls for a balance between 
initial deployment and reserves. In practice, a multitude of 
factors have to be considered before deciding on a partic- 
ular compitment. The approach in this chapter provides us 
with a set of mixed strategies but does not consider intan- 
gJiable factors like world politics, national €conony, 
survival factor and So on. 

Stability has been shown to effect trajectories which in 
turn effect battle outcome. Recall from Chapter IV that 


there are some trajectories which represent speedy and 


complete annihilation of one force; hence it seems reascn- 
able that the side that is tipped to win the battle will 
want to operate onan unstable trajectory. But to what 


extent can one exploit the stability behavior of the system 
to influence battle outcome? Obviously there will be prac- 
tical linitations; an important one of these is total avail- 


able resources. 
Ae PROBLEM STATEMENT AND AFPROACH 


The problem statement is as follows : 


Given total defense resources Q, , Q, for x and y 
respectively, what 1s the optimum set of strategies for 


lnitial force ccamitment, X and Y? 


oa 


We begin by treating this as a 1*1 problem at the gstra- 
tegic Jevel. The dynamics of the problem are thus governed 
Ey eguation 4.1. Both sides are assumed to operate initially 


at equilibrium with constant replenishment rates given by 


ry 
if) 


vein bY 
are aX) (eqn 5.1) 


Y(v + cX) + dX 


Yn 
I) 


Since both sides have limited defense resources Q.0 ye 
the replenishment rates versus time may be as shown in 


Figure 5.1, where Qo. = rT. and Des = st. 


replenishment 





Figure 5.1 Replenishment Versus Time. 


The next step is to select some suitable form of payoff 
function which is to te optimized for a certain choice of X 
and Y. The payoff function (from X to Y) has been chosen to 
be 


x57) @ Lg Ly 


where L,, I, = Total losses for x, y at battle termination 
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As each side runs out of resources at different times, 
the simulation 1S conducted in stages. The total losses are 
determined by simulating the dynamics of the system until 
one of the force levels drops to ten percent of its total 
resources, Q. 

If X and Y are assumed to be chosen from a finite set of 
values, then for é€ach pair (X,Y), one A(X,Y) can be 
obtained. A payoff matrix can be formed and the problem can 
ke treated as a two-person game. Based onthe mirimrax 
theorem, there exists a set of optimal mixed-strategies and 
one convenient way of finding them is through the use of 
linear Progiramming. 

Lt ees perhaps worth-noting that the approach is 
computation- oriented. It has been made feasible by the 
availability of high-speed conputers and efficient software 


for numerical computations. 


Po MUL STAaceE BAgo 


Using the above approach, the entire battle can be 


divided into three stages, namely 


(1) Both r and S are nonzero 
(2) One of the r or S equals zero 


(3) Both r and Ss are zero 


anh 
a 
je 


age i 


This stage will be the period from outbreak Oia 
to the time (T,) when one side runs out of resources. It is 
also possible Wijiat® x 0.19. ory < 0-19, before u is 
reached, in which case the battle is over. In general, this 
period T, can be written mathematically as 


T, = Min {T,, Ty 


a9 


During this stage, the dynamics of the system is 


given by the familiar 1*1 systen 


x = -x(u + ay) - by +r 
: kei a) 
ae - YY (es Gx) dx + S 
Khen this 1¥*] system is integrated, just as in 


Chapter IV, the resulting trajectories behave similiarly. 
However, there is a najor difference. Now, we no longer 
have unlimited defense resources, and this stage will not 
ieee cOrLcGvVer. It implies that, unless Q, or Oy 1s extremely 
large* , trajectories which reflect quick annihilaticn of 
miemot the forces are rare. In general, T, and TY are given 
by 


Q 7 * 


a 
a a 


If one of tke torce levels drops to less than ten 
Pemeent of QO, or Qys the battle is arbitrarily considere 
over and the losses are calculated as in Figure 5.2. The 
finish time (FINTIM) is Simply t, the time when x < 0.10, OL 
y< 0.10,- 


ee 


\ 


Ege 


Since either x or y can run out of reserves first, 
the dynamics of stage 2 are governed by either equation 5.2 


cr 5.4 respectively. 


ae OF may be ver large if x or is backed by a 
superpower on os fully committed wo eouEte Wolab alee rapa acre 
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Figure 5.2 Losses at Stage 1. 


xX = -x(u + ay) - by 
(eqn 5.3) 


-y(v + cx) - dx + 5 


“< 
li 


xX = =X (Uu + ay) = 8b woes: 
; (eqn Sam 
y = -y(v + Cx) —aadx 


Unless the battle ends earlier, this period will 


last for D> whicheis inion, 


T, = Max jT,, Tf ~ ae 
During this period, the trajectory will be different fron 
that in stage 1. This is because when r = 0 or S = 0, cne of 


the hyperbclas is shifted so as to cross the origin and we 
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have different eguilibrium points. The trajectory will now 
be influenced by the new equilibrium point. This is illus- 
trated in Figure 5.3. 

Calculations of the losses are more complex than in 
stage 1 since there are now two cases to deal with i.e. o£ = 


0 or s = 0. The procedure 1S shown in Figure 5.4. 
3. Stage 3 


If the battle enters stage 3 without either x < 
weno, OF YY < 0.10, then the dynamics will be dominated by 
attritions since r =s = 0. Eguation 5.5 is now used for 


integration. 


=x(u + ay) - by 


Se 
t! 


foam 5:5) 
ay (y (tecx) = dx 


“<j e 
i 


Again, the trajectory will have to change because 
now both hyperbolas pass through the origin. This is illus- 
trated in Figure 5.5 where we show how the intersection at 
stage 2 has changed. Losses and FINTIM are calculated in 


accorcance with the procedure in Figure 5.6. 


Ce. MIXEL STRATEGIES 


The range 0to QS for both X and Y can be subdivided 
into m force levels. There are m*m pairs of X and Y and 
corresponding number of payoffs, A(X,Y). We thus have an m*n 
payoff matrix having elements A(x,Y). Figure 5.7 gives a 


pictorial representation of this two-person game. 


restrict the 


>In the actual [froegran to 
pe reflect prac- 


ie wish 
Eange Of xX and Y to interyva 3529) 
ployme 


tical liritations in initial a de n 
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Integrate ieee fi atc 
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x a rit, et oN = X 
ey = s(T, + te) + Yé - y y = Ry 2 5 
eX ,.Y ) = L = L. A(X,Y) = Ly, = Ly. 





FINTIM = t + ies 








FINTIM = t + T 
ee 


Figure 5.4 Losses at Stage 2. 


In the last section, the procedure for computing A (X,Y) 
has been described. A simple program can be written to 
compute each element of the payoff matrix. One such frogram 
is given in Appendix H. 


64 


"¢ ad50Ig Buting UOLYISIOJUYT Irpoqs1odyy puv Sc tC} Oo ne Gc INST 


9 am 
“=e © ame o 
“ro ms eee | 
“~~ , 
“ew » 
—-s , 


, “i 





N a 


391qe3sun a 
9TqQeaAs S 
¢ abeqys -#- 
ONAN Q=Awpon 7 CUS ete 


¢ Sets yeaa. 


65 





integrate eqn 5.95 









A(X,Y) = LL - L 
FINTIM = Max jT\, TY} 
ae 






Figure 5.6 Losses at Stage 3. 


There are a few ways of presenting and interpreting the 
payoff matrix. A normal practice is to present it in tabular 
form and consider only pure strategy. Alternatively, a plot 
of A(X,Y¥) as a function of XK and Y could be obtained. When 
uSing pure strategies, it has been observed that the game 
does not always havea saddle point [Ref. 10] and it would 
be better to use mixed strategies. 

In bixed strategies, x and y may play all their ¢trat- 
egies in accordance with aocertain set of probabilities. 
Although in our situation, x and y can only play once, the 
Same ccncept of mixed strategies is still useful. If we let 
es and gq. be the probabilities by which x and y select their 


ith and jth pure strategies respectively, then 
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FiQUIEC Sey Payoff Matrix. 


Fi 2 P57) 
i 5 


In addition the (1,j)th entry of the payoff matrix be 


denoted by ajj, the probabilities can be represented by the 


Matrix below 


Py eZ RBs, «3 Pm 


| a 


The optimal mixed strategy is based on 
criterion. Mathematically, x and y select p, 


the minimax 


and q which 
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will yield U and V as given equation 5.6 and equation 5.7 


respectively. 
eee | min : a - : 
Ps be 441Pi > DeAgaPiscsss De FamPi eae ae 
i=l i=l i=] 


(eqn 5.7) 
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Appendix I describes how the problem of solving for the 
optimal values of p, and he Cie womolimiitoOulinear progran— 
ming fort. The program given in Appendix H also computes 
this optimup set of solution in addition to obtaining the 
mayort matrix. 

The concept of mixed strategies 1s quite intuitive if a 
game is to be played repeatedly. But Since we are uSing it 
to provide us with an optimum set of probabilities of 
selecting the pure strategies, some interpretation is 
reguired. Although the optimum mixed strategies have been 
obtained, a pure strategy still has to be selected and used. 
However, itis important that the selection process should 
be random® according to the optimized prokabilities 
obtained. 

One Simple but valid statistical procedure [ Ref. 11] to 
select a pure strategy from a set of mixed strategies is to 
first plot the probability distribution function. A randon 
humber generator 1s then used to generate a nunber between 
zero ard one. The corresponding value of the strategy could 


then be selected. This procedure is shown in Figure 5.8. 


6The selection process must be random otherwise the 
opponent can select a strategy to improve his outcome. 
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P(X<x) 


random ([-«<-s=-s=20222— a= 
number 


Strategy x 
selected 


Figure 5.8 Obtaining Pure Strategy from Mixed Strategies. 


D. EXANFLE USING KOREAN WAR DATA 


One cf the main objectives of using actual histormeal 
data in a model is for validation. It is important that the 
results obtained using the model should at least be consis- 
tent with actual events. The Korean War has been chosen 
because there was aclear-cut victor during the initial 
Fhase of the war. We consider the period when only North 
Korea and Republic of Korea (South Korea) were involved. 

Before the entire Simulation can be carried out, the 
actual force strategies, fiyhting ability, weapon State, 
etc, have to be transformed into familiar quantities and 
parameters such as Q,, Qye X, Y, a, b, Cc, d,e--, and Soe 
This transformation, together with some background data on 


the Korean War are given in Appendix J. 


1. Fesults and [iscussions 


= = SSS ea Se 


First we examine the resultant trajectories during 
the three stages of the battle which are shown in Figure 


De oe The simulation uses the X and Y which correspond to 
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the actual initial deployment by both North and South Korea 
respectively. Clearly we see that the victor is x, as it 
was in histcry. The result of the simulation also shows the 
three stages explained in the last section. Note that the 
trajectories for the first and second stages are curtailed 
because both sides run out of war reserves. The implication 
is that in practice, the kind of trajectories leading to 
large and rapid changes in force levels are rather rare. 
However, the effect of instability on battle cutcome 


is borne out by experimenting with the directions of pertur- 


bations. Consider the case in which x (North Korea) fixes 
the initial force andy (South Korea) varying the initial 
force levels around the equililfrium point. Inne tl Tew. 9, 


these perturbed points are denoted by points A to D spanning 
across the boundary separating the domains of attracticn. 
From our understanding of the stability and system dynamics 
€ach perturbation will give rise to different trajectory and 
payoff at the end of the Simulation. Clearly, y will want to 
operate at the perturbed pcints A or B rather C or D since 
the former will result in the trajectory for stage one to be 
in a decreasing x direction. Table II shows the variaticn in 
the payoff as the perturbation point changes. When the 
perturbed points are at A or B, the payoffs to x are less 
then those for points Cor D. Thus we have seen how an 
unstable system can be used to inflict heavier losses on the 
Opponent. The more unstable a system gets, the more signif- 
icant will be the effects of initial perturbation which are 
Manifested by initial victory and element of surprise. Since 
some systems with large aimed-fire coefficients tend to be 
highly unstable, we can expect this effect to tre most 
Pronounced in battles involving high-technology and highly- 
lethal weapons. 

The payoff matrix and optimal probabilities Pp, and 


qd. are shown in Table IIil. The results suggest that the 
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TABLE. 11 
Effect of Different Perturbations on the Passott 
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Figure 5.9 Trajectory for Korean War. 
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North Korears should use large initial deploynent. In the 
actual war, North Korea actually deployed almost all of its 
Pguedueronrce ani within afew days capturel Seoul, the 
Capital of South Korea. The payoff matrix also shows that 
no matter which strategy is chosen by South Xorea, wat anes 
bound to suffer much more losses than North Korea. Again, 
this 1S 1n agreement with history since it iS an accevted 
fact that without US intervention, there would be ro South 
Korea today. 

So far in the example, we have always considered the 
Situation in which the equilibrium point (X,Z%) determines 
the replenishnent rates aS given in edguation 5.1. ees 
interesting to investiyate the effect on the nayoft when the 
initial operating point 1S at some other location cther than 
Pome. | Let the meW initial point be at (X,/¥,) and consider 
Geease where X, 15 kept equal to X and only ty LS varied. 
(X,¥) has been chosen to be (6./7,3.0). In Fijyure 5.10, three 
trajectcries correspcndiny to Y oe oe tees) «(0 age shown 
together with the hyperbolic intersection during stage one. 
Basically, the trajectcries correspond to the three stages 
of situlation as before and x 1s still the victor. However, 
both the payoff (Ly-Lx) and finish time are slightly 
Gifferent from operating at (X,Y). Table IV shows thit y¥ 
inflicts mcre losses on x when opé€rating at Y, above the 
boundary curve rather than at (X,Y), but in doing So y is 
defeated faster. Thus depending on his aission, a commander 
can choose to lengthen the battle or initlict sore casualties 
on his opponent by chcosing a Suitable operating point which 


may te other than an €guillbriua point. 
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Figure 5.10 Initial Operating Points at Non-egqailibrium Points. 


TABLE IV 


Effect of Operating at Non-equilibrium Points 


Yj (Ly -Lx) Finish Time Remarks 
(FINTIM) 
eee A a OE eS Below boundary 
2.435 O. Sale At equilibrium 
Cee) 0.293 Above boundary 
Zee oS Oo Above boundary 
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VI. CONCIUSTONS AND RECOMMENDA Te 


Ae. CONCIUSIONS 


This thesis has covered a number of subjects which are 
based on the generalized Lanchester Model. The first part of 
the results has to do with finding the equilibrium points in 
the N¥N system. The Continuation Methods have been found to 
be suitable for this purpose. The advantages of the 
Continuation Methods over numerical technigues are numerous 
and important to our understanding of the non-linear set of 
equations. The method finds all the eguilibrium soluticns 
accurately and does not need good initial guesses. An 
example to compute tke equilibrium solutions of a 2%*2 systen 
is presented along with a way to treat singularity fproblen. 

The derivations and interpretations of the relaticns 
between stability and system parameters form the next major 
portion of the thesis. By considering the simpler 1*1 
problem, a few interesting conclusions have been reached, 


nanely 


(1) Poth the system asymptotes and equilibrium points 
are intrinsic to a system in equilibrium. The locaticns 
ef the equilibrium points on the x-y and €,, €, planes 
ccmpletely characterize their Stabamiaries- CS and Ey 
are the differences in the system asymptotes; 

(2) The dynamics of a system are characterized by the 
phase trajectories which represent the ways a battle 
progresses. Besides stability, the domains of attrac- 
tion also influence the trajectories. The boundary 
curves which separate these domains can be ascertained 


by graphical or backward integration. 


Ts 


The last portion of the thesis integrates the concept of 
equilibrium stability and system dynamics. It relates these 
theoretical concept to operational problems. Two operational 
issues are addressed namely, (1) the effect of varying X and 
Y, the initial force deployment on battle outcome, (2) the 
exploitation of stability to influence battle outcome. A 
methodolcgy which comkines multistage battle simulation with 


two-person game has been emfloyed and the conclusions are 


(1) Initial force deployment, X and Y can be optinrized 
by finding a set of mixed strategies. A suitatle pure 
Strategy can then be selected from the mixed strat- 
egies; 

(2) Instability can and should be used to share the 
course of battle and its outcome. This is particularly 
true in highly unstable warfare which is normally asso- 
ciated with large aimed-fire attritions. Unless defense 
resources are extremely large, it is not possible to 
completely reverse the outcome of a lopsided-battle 


where cne side is much stronger than the other. 


As far as military commanders are concerned, the above 
conclusions suggest two things. Firstly, depending cn the 
relative strengths, Lt As not necessarily true that 
deploying the largest possible force will bring victory, 
reduce loss or even buy time. There is an optimum way of 
depioying available forces. Secondly, if a war involves 
large aimed-fire attritions due to weapons like aircraft, 
Missiles, tanks, artillery, naval bombardment, etc., then 
initial victory which could perhaps be achieved through a 
preemptive strike certainly affects battle outcone. 


Ts 


Beh RECOMMERDATIONS 


1. Iransformation of N#N Eroblem into the #1 Probier 


It has been mentioned that the sinplicity of the 1*1 
problem can be attributed to the Simpler mathematics 
involved and our ability to draw and visualize two- 
dimensional pictures. Despite its Simplicity, it does share 
many of the properties with the N¥N problem. Considering 
these factors, it seems logical that an attempt should be 
made to find the 1*1 equivalent to the N*N system. Another 
reascn is that there is much to be gained in terms of 
Savings in computational effort by going to thera 
equivalent. 

Of course, the "equivalent" system will not be 
expected to be identical to the N¥*N problem in every asfect. 
One can only hope that it 1s equivalent in some sense, for 


example 


(a) Preservation of stability characteristics and 
dynamics; 


(b) Preservation of mixed strategies. 


One way Of transforming the N*N system into the 
equivalent 1*1 system is to equate losses in both systeis. 
The eguivalent system parameters are obtained by using rela- 


tions such as 


Manca = 202 pag 
1 J 


ly Ys = coordinates of the equilibrium point in 
the N*N problem 
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eq 41 
i 

‘eq 7 DY; 
J 


The results of the preliminary studies suggest that 
this method of transformation can preserve some stability 
characteristics. The possibility of using the eguivalent 
system to obtain the mixed strategies should not be 


dismissed until further studies have been conducted. 


Ze ‘Time Variable Replenishment Coefficients 


The replenishment rates used in the thesis have 
always been assumed to be constant. In actual wars, constant 
replenishment rates may not ke used by either side; at 
times, it may not even be possible to do so. It would there- 
fore be interesting to study the cases which involve tine- 
varying replenishment rates L(t). The choice of r(t), for 
example periodic, non-periodic, ramp, etc. will depend on 
how well it represents practical replenishment rates. 
Whether a mathematical tool can be found to cope with the 


added complexity also has to be considered. 
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APPENDIX B 


NUMBER OF EQUILIBRIUM SOLUTIONS OF THE N¥N PROBLEM 


We first consider the 2*2 and 3*3 problem and extend the 
results to the N*N problem. 


1. 2*2 Problem 


In general, each trivial solution leads to a unigue 
equilibrium solution in a continuation process and there 
will be egual number of equilibrium and trivial solutions. 
Secticn 4 discusses what happens when there is degeneracy. 


The trivial solutions are obtained by solving 


Z, (uy +4, 123441274) 


= 0 
Z 2 (uy*ay123%82774) (eqn B.1) 


mgs 3 171°°2172? 


eT.” 


At first glance, there would seem to bea 27 trivial 
solutions corresponding to the number of ways of making the 
lefthand sides of equation B.1 zero. Each lefthand sides can 


ke made zero by either making Zz; or the terms in parenthesis 


equal to zero. But closer examination reveals only six 
allowed cases, in non-degenerate cases, corresponding to (1 
+ 2° + 1) = 6 solutions. Table V shows how these cases 
arise. 
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TABLE V 


Trivial Soli@wien for 2° 2 8emeeen 





Case Zz. which ane eieGiem as Remarks 

i Z4, 25, 23, %y 

2 Zo 23 

: rials 2° = 4 cases 

4 Za» Zo 

5 Zn, 2y 

6 —_ all terms in parenthepgeme 

equal to 0 
7 Ze degenerate case; 


see section 4. 


2. 3%3 problem 


Here, the trivial sclutions are obtained from 


oy 21174 eae ao se 
Zo (U5*ao424t8r92c48526) = O 
oz (uztazs2ytazoZ%etaz226) = 0 

- 0 


Z gC g ty 2476242943423) 
Zp (Uic+C1 72" Gap) <7 


Ze (UG*C) 521 "Cee ees) 


Again, there are some cases which are not allowed 
because cf inconsistencies in non-degenerate cases. Table VI 
shows the different cases. 


In this case, the total number of allowed cases iS 


2 BNR 
p+) +G) +2. 20 
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WD lala ol 


Tio me moet clone O37) 3)cProblem 





Case cs which are made zero Remarks 
ik Zyy Zaeeey Xe 

2 Zyy Zy 

: ea: 

: m2) “6 () = 9 cases 
5 Zny 2y L 

6 254 Z¢ 

i az 

8 Za5 2, 

, #39 °S 

10 a3) Ze 

a me eee ey 7s 

EC an eee > EG 

i. ome > 5) oe 

Pe “2 _ 2 - 

17 - a 3 _ 

a 

r eee 3 EMS 

a mea? >, 

20 es all items in parenthesis 

equals to Q. 
21 Cea 2G degenerate case; 
see section 4. 
BE zy degenerate case; 


See section 4. 
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3. Extension to N*N Problen 


Now we know the type of allowed cases, it is easy to 


recognise the pattern of results. The pattern looks like a 
Pascal triangle with each element squared. 


Number of 
N*N eee . 
Equilibria (N..) 
a 1¢ at = 2 
#2 1¢ + 26 5 a =e 
3%3 12 «ge See = 20 
2 i 

4*4 p24 gt + 62 + 4% 4 1 a 

2 2 2 2 2, 34 = 252 
5*S 1° + So te ee 


For N¥N in general, 


ue Can also be written as 


i=0 


can be found in [Ref. 12]. 
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For the 2*2 problem, a degeneracy can arise when Zy 
(or any cther Z 5) is made egual to zero. The other three 
lefthand sides in equation B.1 are made equal to zere making 
the terms in the parentheses equal to zero. 

In general, such a case does not correspond to a new 


trivial solution because there 1S an inconsistency when z., = 


2 
“u,/c., and Z, = ~U,/C,5- However, it seems that Fy making 
U/C, = U/C 556 the inconsistency no lonyjer exists and 


there will be an infirite number of trivial solution as long 
May ee 3 To. |CT 
turns out that the nrumber of egquilibrium solutions still 


aS Z,, 2, are chosen to satisfy u + a 
remains at a maximum of six (disregarding the case with 
infinite number of equilibria). Similarly for the 3%*3 case, 
even if there is degeneracy, the number of equilibria will 
net exceed 20. 

That the above is true can be shown by considering 
the 1%*1 problem. In this case, the two non-degenerate cases 
correspond to (a) x = y =0 (b) v + cx = 0 and u + ay = QO A 
degeneracy can occur if v = 0 or u = 0 in which case there 
seems to be an infinite number of trivial solutions lyirg on 
the y or xX axis respectively and hence an infinite number of 
equilibria. But when the actual hyperbolas are plotted, 
there are only two intersections and hence two equilifrium 
points. Furthermore, when the Continuation method is used to 
find the eguilibria, there are only two equilibria irrespec- 
tive of whether the trivial solutions are chosen to be 


degenerate or not. 
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THE SINGULAR CASE IF REFLAG OR ERFLAC 
STEP: 


1 
CALL TREAT (CREFLAG,ERFLAG,CHFLAGI 


ME AK MA ME Rh KK SSS SAKE KKH 


WHILE T<il OR CORRECT CR STEP OCES NOT RETURN WITH 
TREAT 
[oe 


ERRUR FLAC IER 
RMR BG MK BRE RGB RK RRR RK RE EK KK EH 
SET UP_THE SET OF DIFFERENTIAL EQUATION ANDO 


CALL LEQT1F TO SOLVE THE SYSTEM AX 
$I RR BRR LR RR RR SK 


Leek FRE SHE DED ER ARE RS CHEE ERE REE SERRE ERS REEES EEREEE 
PREDICTOR 


YVOVOWN VUOOYO VYOVUVQVUOY 


~a .J a) —J a. 
CL kel hy Uy 1 
OA WA a 
we wet te = 
om “Qe — — <I 
PIN Nw ~ ~ ~—ac 
et ret et ee >< x KAA 
mh >E OE OMS He + t+ +x 


wet et eet ee I {Ian w=" em OT" Me 
DW + OHNO +ANAAFONNO +a 
trretro Eee eo a 6a oe eee 
ae ee EHO) HHONNNANQOANOTU, ow 
HAN) See wer ese ee er er er ee Ow 4 10. 
wwe wr es IO AT ODA TOOQDOVVO Dre i} @ 
POCO ac at iy Wn ee Ned gee ep ee 
Od et a a ee ag Ne ey em a RN me ER mE BAL) | ==) 


WU BEANO RAND PANO PANO INS 


—$ | —E-—m=-m & & © SF S&S SF & F RS & & SF & & B&B C~wwe4lL 
ANID SF tt AMIN AUN NID MO OP PP dt OS Fe 
ad et et eget Set ee? ee ee See ee OE gy ened eae eee eg ene he? eet “et al eee 
ONQAQAagaaaaadaadgaaadaanaowaag Orr 
mMmMUmMmuaqidqad<qadqaqtataatttyw Ue 
= 

CO CQ 
© YO 


LOZ 
101 
C 


i) 


CK ee RS RR eR a a a oe a a a a a 
CORRECTCR STEF: 


C 





CR YR SS A a HK BK KKK 


CALL ZSCNT (FCN»NSIGeNy 1TMAXe Ty XE NCRMy HK, LERD 
OMAX=25C.O0DO#DMIN1(DX(12,DX¢2) 9 CXC 2) 9DX(4)3 


C 
C COMPUTE DIFFERENCE BETWEEN PREVIGUS AND CURRENT X 

DO 104 ks 

bck =D ABSUXT(K)-X(K)) 
104 CONTINUE 
C CALL SUBROUTINE TO DETECT FAST CHANGING COMPONENT 
; CALL DETECT(DX »,REFLAG, ERFLAG) 
Cesssessecsseersoanete sate tess avaedan needs sae sonst texe es 
C MPUTE ERROR IF T 1S NEAR 1.0 
ee ee ge ss 6 444 oe doe Ko ES Ha R HEHE 
C 
~ Fl==$X(1)*(U (01) tA (Ll) L)*X (3) 4A (1, Z)*®X04)) 

& tee eX (3IEB CL, 2)4X04) 
F2==-X(2)%(U (2) A (2) 1)*®X (3044 4252)*X(4)) 

& +h Come 1) *X (3 )=B 02,2 78X04) 
F3=-X(3)#(U 632 4C (1, 1)X 4104C (2, 1)*X(2)) 

g meee 1,1 )*Xi1)—0 (2,1) #XC2) 
oR Ca)ob Lod eK CLO Card FREE d 2)*X(2)) 

& Hae (1.2 14X(1)-D(2,2)%xX12) 
ERR=CSQRT (F 1#*¥ 246 24K 246 3KS24F4% 42) 
WRITE(L »997 1T(1),ERR 

1031 CONTINUE 
RR MR A MK RH KSA HH 
QUIP UT VALUE CF X VECTOR WHEN T IS NEAR 1. 
CR RRR SKS RAG RR VKH X MK RE Me Me ee eR eR RH KK 
é 
LF ( .NGT o(T (1) sGE« (0.9900)) 160 TC 865 
WRITE(1 »955 IT(1) 
WRITE(1 ,-999 34K pK (CK) »K=1 NI) 
865 CONTINUE 
CK RR BK Me Re DK Ke HK EH 
C CHANGE DOUBLE TU SINGLE PRECISION FOR PLCITING 
CM A HE ER BK KSSH HVE EK SEM 
E 
TA(CTRI=SNGL(T (12) 
XA (CTR) =SNGL(X O12) 
XB (CTRI=SNGL(X (2) 2 
XC (CTRIESNGLO(X (3) ) 
XD (CTR) =SNGL(X (4) 3 
COUNT=CCUNT+#1 
CTRK=CTR#l 
GQ 10 5 
: CONTINUE 
CR RS SKK RK KR RM RK RK RK SK RE 
C COMMENCE FLOTTING ROUTINE 
Ce eA Be Me Mee A HH Me a Re Me Re HH Re ee eH HK HK HK HK KH 
C 
CALL TEK618 
c CALL COMPRS 
CALL NOBROR 
CALL BLOWUP(0.4) 
CALL AFEAZD(12.0, 80) 
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APPENDIX D 


DERIVATIONS OF THE RELATIONS BETWEEN aa 7 AND STABILITY 


1. NewWtead Stabila, 


In this case, _— at = 0 and equation 4.5 reduces 


xe, 

Ny (X-x,4) = (XY-X 01% e7) e Laiaaee sar 
It follows that both sets of hyperbolas merge inte one and 
all the points on the commen hyperbola are equilibrium 


points. The lefthand side of the second condition in egua- 


tion 4.4 can be manipulated as follows 


(utay J (vtcx,) - (b+ax,) (dtcy,) 


i 


ac [(nyt¥g) (H2*X,) - (Hy txXQ) (ngty,)| 


ac [x,(ny-np) * yg(uy7¥y) * nyyy - n2¥1| 
= We 


This result implies that the constant term of th2 character- 
istics polynomial D(s), is zero ; hence one of the eigenva- 
iues, S_ €quals (Comer. Factoring out the characteriaue 
polynomial, we get the other eigenvalue as 


Oe mals | (uray,) ~ (veex,) | 
(eqn Daw 


- [acnyty,) + c(u9*%Q)] 


i 





FOL Ppormes On the first gquadmpant hyperbola, x, > salle and ie 


> Sais os therefore s, <0 and neutral stability exists. 
Conversely, S., > 0 on the third quadrant hyperbola which is 


therefore unstable. The resultS are summarised as fcollcws : 


(1) We Se = Sn , infinitely-many equilibria 
€xist as points cn the two hyperbolas on which x = y = 
C; 

(2) The first guadrant hyperbola is neutrally stable; 
(3) The third guadrant hyperbola is unstable. 


2. Intersections in First and Third Quadrant 


The proofs for the following results are given in 


mes Section : 


(1) when both equilibrium points are on tne first 
quadrant hyperbolas, one is stable and the other 
unstable 

(2) when one eguilibrium point is on the first gquad- 


rant hyperbola and the other on tne third, both can be 


The straight lines given by equation 4.7 are plotted 
Om the c,, Ey plane as shown in Figure D.1. It also shews 
the corresponding regions on the ¢€,, ¢€ 


plane as x, and Te 


MC 2 


vary. 


Let (X46 Yo.) be the first eyguilibrium point on the 
first quadrant hyperbola. The first stability criterion in 


equation 4.4 1s autozatically satisfied since 
(utay,) + (vtex,) 
7 a(ns+y,) + C(u5+*x,) 
> 0 


Ou x9 > “U5 and Yo? > “Ny 


Ke 





Figure D.1 Effect of Varying X and Y., ME wee Plane. 
x i 


2 


Consequently, the stability of (Xa eee is solely deter- 


Mined by the second condition in eyquation 4.4. The second 


condition can be rewritten in the following manner 


Cutay)(vtcx,) ~ (btax)Cdtcy,) 


aS | x(n yon) + Yoho =e gang * ny (ug-Hy) | 


ac “E(x, tu? + Eni ty,) | 


a) 
S 
or Ey Sex (M1 *Y,) (eqn Dag 
. (u, +x.) 
1 “e 
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Fox (X 16 Yeas equation D.2 represents the boundary 


line of stability and this line has a slope between ny (%, 





1 
+ Wj) and (Yo, * "7, - Tt is shown in Figure D.2. 
Pagure <Dii2 Boundary line of Stability on Ee = Flane. 
To investigate the stability of (X55 oda we 


substitute them into equation D.2 and obtain the condition 
for stability as 


u Ge a) 
: > —— (ean D.3) 
2 Cu *X 67) 


100 


Equations D.2 and D.3 are only “Gietecient in eee 
directions c£ inequality. This means thaw on one side of tae 
line Ca Ee (,* Y, 7 0ca ee ee one equilibrium foint is 
stable and the other is unstable. On the other side of the 
line, the opposite ccnditions exists. 

Note that (X.5,, Y,5) May or may not satisfy the 
first condition of e@wation 4.4. This Condepesoneins 


a(n, *Yo) + C (ute +X, 9) Sb (eqn Doge 


From Figure D.1, eguation D.4 and earlier results the 


follewing deductions can be made 


(1) In the first quadrant of the Cae, plane, ea OF 


xX 5 > aloe yo 2 “3 hence equation D.4 is satisfied. 


, WY ; , 
The region labelled YH); in Figure D.3 1S stablewiae 


(X58 Y.>) but unstable fer (X.a¢ Yo ,)3 


(2) In the third guadrant of the Ele Ey plane fut 
ketween the two lines where x > 0 sand Y., 2a 


e2 I : 
eguation D.4 is again satisfied. The region: labelied 
AA is also stable for (X50 1) but unstable for 
(Xue Yo) 9 
(3) In the second yuadrant of the owe, plane, %X.,< 
— ie ey < -n. and e« < 0, equation D.4 iS not Sages 
L e2 Ie X a 


freien Ss OmeO ti (X10 Yo.) and (X.5- py, are unstable. 
The region is labelled Mm Pagure pss: 
(4) In the fourth quadrant of the Ee - plane, equa- 


tion D.3 is not satisfied ; so (x is unstable 


ees ey 
but (X.7- Ya) 15 s@apie- This region is labelled 


in Figure D.3; 


eueMmeecaanete Se eecaeh-4 
tA eed ede ten! rr ho 


(5) In the region labelled BeHBORe, eguation D.3 is 


ww sy | ttt wt LL 
“Tr rr tits fF PP ee 






Nera ete eertenhon 


not satisfied ; so (X,5- ¥,3) is unstable Fut "(ge 


yoy is stable. 


(Xe) Yeo 
WA. +s unstable; 





' (Xe2> Ye2) 
Stable 
Non as above 
=o both 
/ = aie unstable 
o AEE Sth 
TT 
a eo FeeEmaaenmert | ic Cerys) 
UCC OSaaneSeeanee. eligean €k 





oa aq 
oS tale) ell - ‘ ) 
| aml eo 2 
a unstable 
: wagoocasn 7 
: Safe aaa HH §=6aS above 
Se Ee aes f + 
eee, 
SS 2 coe) eels (n,+y_,) 
SS PP elmitejaiere| ae | Be oe eae. 
Pammee7Z Nala: i.) 
a aaaalcied eet es 
S feet * HS esa 


Figure D.3 Regions in €,, a Plane. 


By carefully noting the signs of (X,5, Y,5) and the 
various regions in Figure D.3, all the previous deductions 
can be ccmbined ; we conclude that : 

(1) When both equilibrium points are on the first 
guadrant hyperbola, one is stable and the cther 
unstable; 

(2) When one equilibrium point is on the first guad- 


rant hyperbola and the other on the third, both can be 


unstable or one will be stable and the other unstable. 
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3. Iwo Equitjibeaa in the Third OUa@iaame 


The simplest way to show that both eguilitrium 
points in the third quadrant are unstable is to refer to 
Figure D.4. Clearly, we must have x < - vfc and y < - w/a 
for koth equilibrium points. In that case, the first 


Stability criterion in equation 4.4 is not satisfied since 
+ + 
(u ay) (v+cx ,) < 0 


Therefore, both equilibria are unstable. 


~ x owe 
$e 
_ — — ——. -_ - = —_—- «+ - oe eee ae = = «@ ee ea 
- = 2 «@ A ibe ee 








Figure D.4 Two equilibria in the Third Quadrant. 


4. Repeated Equilibria 





This case corresponds to operating exactly on the ¢« 


MY 
> : i 
c + nj7(x uy) on the €,, Ey plane. The procf is 


obtained by substituting Xoo = Xo. © Yoo = Yay into equation 


4.6 and solving for o 
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= = ee : 
mez Jelw C, (xei 7H) 7 7y 
; (eqn D.5) 
a es : 
elm seul E (Ye, *y) Mi 
Rewriting equaticn D.5, we have 
Ec rees ~y el Tew fox 
cies x el Wx © N1*y 
subtracting one from the other, 
= 2 é 
ce tenl orl . elie, cae 
Ceol 
6 fe sel ol" ec ion) 


y x 
(X41 7H)) 


Thus we have shown that the case of repeated equi- 
libria corresponds to points on the straight line indicated 
pieragure D.3. 

Eguation D.6 can be also obtained by setting the 
second stability criterion equal to zero ( sée eguation D.2 
)- That 1s equilavent to saying one of the eigenvalues, s, 
is equal to zero. the sign of S, is determinei by censid- 
€ring the first stability criterion. Following the same 
argument as in neutrally stable case, we can prove that 
repeated equilibria on the first quadrant hyperfclas are 
neutrally stable. On the other hand, repeated equilibria on 
the third quadrant hyperbolas are unstable. 
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APPENDIX E 
RELATICN BETWEEN STABILITY AND Ee Ey PLANE: VERIFICATIONS 


Scme representative points on the Ee 2 plane are 
chosen and the corresponding stablities of the equilibria 
calculated. The points chosen are marked F, G, H, M, N, P, 
Q, T and W in Figure E£.1. 


Cy / f 
/ ra 
H a 
: / N 
p [ a ee 
—— o — 
IG.” oan 
F / - Gr 
——_— $f >” 
* a a / 
‘ a i 
i / 
a 
a [ 
on Ve 
W / 
/ 
, eT -Q 
/ 
/ 


Figure E.1 Experimental Verifications. 


1. FErocedure 


Having chosen the points on 
have to work backward to obtain a, 
Suitable (X17 Yo ,) are then chosen, 
tron CL I, S, X 5 and a (X 5 
directly from equation 4.6. From all 
eigenvalues can be calculated and 


theory. 


¢- Results 


The results are tabulated in 


with tke theoretical results given in 
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Cores  y  Piahe, we 


b, Cac u and Vv. 


followed by a calcula- 


Y 52 
these parameters, the 


) can be cbhtained 


results compared with 


Table VII. They agree 


Figure £E. 1. 
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Mees DIX I 
TRANSFORMING MIXED STRATEGY PROBLEM INTO LINEAR PROGRAMMING 
iiewoayvort tinetron has been defined to be 
ewer s Li = Lig 


Meese ects fis ODtimum mixed strategaes which veild 


mM 
“}> ; 2 i 
V = 13° je a>: ieee 4m 3 45 


J! 
where 
i = payoff to x when x adopts ith strategy and y adopts 
jth strategy 
we Secu dmrevmcipate y Selects jth strategy 


Siegect to 


iiiwmeche Original problem becomes 


minimize Vo 


ae 


Subp ]eCcu sae 


Dividing the constants by Vo (>0), we have 


ip Sa 
beeen 


Ne | 
Sees 
Je 
1 m 
Foes 
=] 
q . 
ket) Ca _J and since 
J 0 
min Vo = max : = max LQ, Q, Sey esse oe Ql» 


the problem can be written as 
Max Qh z= [Q, 4 Q, eee Qn 


lus 


Swine CU mee 


meeetesolvyine the LP problem, the optimum strategies for y 
* * 
is given by a5 = Q; Vo Some constants, K = min(a, 


have been added to =e to ensure Vy > Ao ethics: 1se-done, 


)| could 


K has to be subtracted from the optimum value obtained by the 


ieee that 1s 


where VOT eenomvalue Of the Gamer 
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APP ERD I Xx 


BACKGROUND DATA ON KOREAN WAR 


PERIOD CONSIDERED ™ * 235 June 50> to a ae 
(No American ground involvement 
as yet) 


TYPE OF ENGAGED 2) ieee Predominantly land combat 
GENERAL STATE OF REFIE ESS 


NORTH KOREA: Well prepared by 1950; arms 
build-up and trainingees 
troops since 1945; many 
military leaders and combat 
personnel were war Vete vam 
PMehtineg gm Chime 


REPUBLIC OF KOREA: By 1950; a small defense 
“~ force began to take shape 
through American aid; 


training only started arogmm 
1948. 


SOURCE OF DATA : a) Appleman, R.E., United States 


“Army in the Korean War, Department 
Of the inmy sanso! . 


b) Montross, Lynn, U.S. Marine Opera- 
tions in Korea, U.S. Marine Gomes 
sam . 
RELATIVE SiRENGINAD : 
NORTH KOREA REPUBLIC OF KOR 


Total strength = Q. = 135,000 men Q, = 95,000 men 


Tanicss:: LS0 eran 

Mec IIe jeeces = Lowe 700 

Alret@amt 

(2) fighters - 40 no combat aircraie 
(11) attaek  bembeis a (22 trainer, 4 

(111) reconnaissance - 10 auxiliary; no pi lGmm@ 


OF Cehieeor NDING PARAMETERS USED IN MODEL 


a = 0.7 c= fC 
mee 0.4 d = 0.6 
u = 0.15 Ver (a2 


euin © = >. 5 0 ene m 
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